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We address the computation of physical observables in graphene in the presence of Coulomb 
interactions of density-density type modeled with a static Coulomb potential within a quantum 
field theory perturbative renormalization scheme. We show that all the divergences encountered in 
the physical quantities are associated to the one loop electron self-energy and can be determined 
without ambiguities by a proper renormalization of the Fermi velocity. The renormalization of the 
photon polarization to second order in perturbation theory - a quantity directly related to the optical 
conductivity - is given as an example. 

I. INTRODUCTION 

The role of many body corrections to the physics of graphene is at this point uncertain. While in the first transport 
I 1 experiments electron-electron interactions seemed not to play a major role [3-0] more recent measurements 0-0] 
and the observation of the fractional Quantum Hall effect [l0|, LLL| indicate the possible importance of the Coulomb 
, interaction to the intrinsic properties of graphene. 

The problem of the nature of the interacting system was addressed in detail in the early works on graphene both 
in the weak [l2l - fl5| and in the strong coupling limit [l(l-[l8j . The main issue was to determine the infrared nature 
of the system -its Fermi liquid character and its associated physical properties-. In the weak coupling regime [l2[ 
with a perturbative renormalization group (RG) analysis it was shown that intrinsic graphene behaves as a strange 
Fermi liquid in the sense that all possible interactions are marginally irrelevant but the inverse electron lifetime grows 
linearly with the energy instead of quadratically [131 ] . An upward renormalization of the Fermi velocity at low energies 
was also predicted what in turn implied a downward renormalization of the effective Coulomb interaction (g = 
to the infrared. Experimental indications of both the linear inverse lifetime of the electron 0, 0] and of the Fermi 
velocity dependence with the energy [9( have been recently reported. The results of these early works have been 
reproduced and pushed forward under several approaches after the synthesis of graphene [l9i - [23j . 

A very important aspect of renormalization was left aside in these previous works: the renormalization of the 
parameters of the theory that allow to determine the physical observables. We will here follow the standard quantum 
. field theory (QFT) approach to renormalization aiming to give a prescription to calculate unambiguous, cutoff- 
independent observable results to any order in perturbation theory. We show that a proper renormalization of the 
Fermi velocity renders all physical quantities finite and unambiguous to second order in perturbation theory. We 
work out as an example the renormalization of the photon self- ener gy, a q uantity directly related to the Coulomb 
interaction corrections to the optical properties of the system [2ll . l22l l24 - |2q ] . 



II. RENORMALIZATION IN QUANTUM FIELD THEORY 

Ultraviolet divergences arise in QFT due to the singular behavior of the fields at very short distances in real 
space - or at very large energies in Fourier space-. Renormalization [53, [H| is a prescription to get rid of ultraviolet 
divergences and construct sensible models where physical quantities can be accurately computed. The idea is that 
the ultraviolet divergences can be canceled by adding counter terms to the Lagrangian what amounts to a redefinition 
of the parameters (mass, coupling constant, wave function) of the theory. The process is usually done order by order 
in perturbation theory. If done appropriately at the end one finds finite results independent of the regularization 
procedure in the computation of physical observables. The prize to pay in this process is the necessity to fix the values 
of the renormalized parameters at a given energy from some well-chosen experimental inputs. This "renormalization 
prescription" is crucial in the process and the basis of the later renormalization group (RG) developments in QFT. 

The deduction of the experimental values of the parameters from a given experimental measure often involves 
phcnomenological assumptions. Very good examples of these difficulties are provided by the interplay between the 
accurate determination of the fine structure constant of QED - that often is done from solid state measurements of 
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FIG. 1: Tree level Feynman diagrams (see text for details): (a) Electron propagator (b) Photon propagator (c) Interaction 
vertex. 



the electron precession in magnetic fields - and its feedback to determine the anomalous magnetic moment of the 
electron with the actual precision of better than one part in a trillion [29| . 

In the condensed matter applications the difficulty increases due to the fact that we do not have "scattering" 
experiments involving the asymptotic states of the fields but transport measurements that are influenced by all kinds 
of extrinsic factors (disorder, doping, substrate). Yet the renormalization program can be adapted to condensed matter 
systems that admit a continuum effective description, graphene been one of the best examples. The independence of 
the observable quantities on the cutoff guaranteed by the procedure makes it irrelevant whether or not the cutoff has 
to be taken to infinity or to a finite value defined at high energies as the inverse lattice spacing. It ensures that the 
observables of the effective low-energy theory do depend on the high energy only through the renormalization of the 
effective parameters and do not have an explicit dependence on high energy quantities. 

Next we will see these words at work in the concrete case of graphene physics. The program consists of identifying 
the primitively divergent Feynman graphs, adding counterterms to subtract the divergences, redefine the parameters 
of the model to absorb the infinities, and fix the finite parts of the vertex correlation functions by a renormalization 
prescription. For this last step we need a number of external conditions (observable data) -the renormalization 
conditions- equal to the number of parameters to be renormalized. The original RG equations in the QFT approach 
were established to demonstrate the independence of the observable quantities on the renormalization prescription. 
Any two set of experimental data will give rise to the same result for a cross section. 

III. RENORMALIZATION OF THE GRAPHENE MODEL 

A. Definition of the model 

It is known that the low energy excitations of a graphene sheet around a Fermi point are well described by the 
massless Dirac equation in two spacial dimensions [30, |31| ■ The density of states at the Fermi level vanishes and the 
Coulomb interactions are unscreened [l2j ]. The non- interacting model considering a single Fermi point is described 
by the Hamiltonian 



H = hv F J ^rVKrb^VKr) , (1) 

where i = 1,2, 4>(r) = ip + (r)j°, and the gamma matrices can be chosen as j x = <J2,ly — = 03 • are the 

Pauli matrices and vf is the Fermi velocity (to be defined unambiguosly in the next section). 

We will follow ref. [lij and model the electron-electron interaction as a density-density interaction mediated by a 
scalar potential. The (instantaneous) Coulomb interaction can be described by the scalar component of the gauge 
field: 

H in t = e / d 2 r^(r) 7 V(r)A (r) , (2) 



with our choice of gamma matrices. In order to define the renormalized theory we need a Lagrangian, a renormalization 
scheme (i.e. a regularization method and a set of renormalization conditions), and the experimentally measured 
parameters associated with these conditions. 
The Lagrangian is 

C = I d 3 kip [7% + W7 • k] V - e^7>A). (3) 



The Lagrangian ([3]) contains four quantities that can be redefined: the velocity i>, the parameter e in the interaction, 
and the electron and gauge field wave functions. 



b) 





FIG. 2: Primitively divergent Feynman graphs in QED(4). 



FIG. 3: Tree level counterterm associated to the electron self-energy. See text for details. 



The electron and photon propagators in momentum space Fig. Q] (a) and (b) are given by 

7°A;o + v-f ■ k 



Go(fc°,k) =i- 



-(feO) 



w 2 k 2 



(4) 



(5) 



The tree level interaction vertex (Fig. [T] (c)) is F° = — icy . Notice that the gauge field propagator does not have the 
canonical QFT dependence (1/fc 2 ). The reason is that, although the electrons are confined to the two dimensional 
plane, the electromagnetic field lives in three spacial dimensions. The bare propagator in {5J makes the graphene 
model different from QED(2+1). In particular the interaction term in the Lagrangian is scale invariant: the coupling 
constant of the present model is dimensionless unlike that of QED(2+1). 

The model was shown to be gauge invariant and renormalizable and in [l2j . The renormalization functions can be 
defined from the self-energy and the vertex as : 



Gq 1 - E(*°,k) = Z7 1/2 (fc°,k)[A ; 7 - Z v (k°,k)vj ■ k], 



(6) 



r = Z e e 7 °, (7) 

It is important to note that when renormalizing at a given order in perturbation theory, the vertex functions (am- 
putated one particle irreducible Green's functions) directly related to the observable quantities have to be computed 
as a sum of all the corrections and counter-terms up to this order. It does not make sense to renormalizc a single 
diagram. 



B. Renormalization of the theory to first order in the interaction 

In QED(3+1) there are three primitive divergent diagrams shown in Fig. O As there are three free parameters in 
the model: the electron and photon wave functions and the coupling constant, the model is strictly renormalizable 
and when computing a higher loop diagram we get higher powers of the logarithmic divergence. 

In the model of instantaneous Coulomb interaction for graphene the only primitively divergent graph at the one 
loop level is the one corresponding to the electron self-energy in Fig [2] (a) and of this, the divergence only affects the 
spacial part of the momentum. The result of the computation of the diagram with a hard cutoff is: 

E£> (k) = - |trr • k (- log ^ + 4 log 2) . (8) 

The electron self-energy can be made finite at this order in perturbation theory by including at tree level a counter- 
term of the form depicted in Fig. [3] with the associated Feynman rule: 

O k ) = § "7 • k (log A 2 + 4 log 2 + F A ) , (9) 
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As expected in a renormalizable theory, the counterterm has the same operator dependence as a term in the original 
Lagrangian. Since the only requirement to impose on it is to cancel the divergent part of the given diagram it 
contains a momentum-independent arbitrary finite part F\ to be fixed by the renormalization condition, i. e. by 
an experimental measure that allows to extract the value of the two point function at a given momentum k R . This 
condition introduces the scale that settles the apparent dimensional mismatch in (|9]). 

Summing up the contributions of the tree level plus the two Feynman graphs of order g (Fig. [2ja) and Fig. [3} the 
two-point function is: 



Gl\k°,-k.) 



(7% 



- v*y 



l-^(logk 2 + F A ) 



that can be written as 



with 



G- v 1 (k°,k) = -t ( 7 °fc + t;(k)7.k) 



«(k) 



1 



|(logk 2 



-Fa) 



(10) 



(11) 



(12) 



The two point function (fTT|) has the same form as the free one with a k-dependent and arbitrary parameter v(k). 
The last step of the renormalization program is to fix the arbitrariness with a renormalization condition. For this we 
need an experimental measure of the Fermi velocity at a given value of the momentum k R . 

The experimental determination of the Fermi velocity of graphene [32| is an important -and elusive- issue similar to 
the determination of the fine structure constant in QED(3+1) and we will discuss this issue further in the discussion 
section. As in the precisson tests of QED, each comparison between theory and experiment can be seen as an 
independent determination of vp. To exemplify how the renormalization works we can take as an example the 
experimental value of the Fermi velocity given in [331 ] : 



v(125meV) = 1.093 • 10 6 m/s = v F . 



(13) 



With this condition we fix the value of choosing the bare velocity to be v = Vp. The physical Fermi velocity will 
depend on the energy at which it is measured and on the renormalization point kp> (125 meV in this case): 



v R (k) = v F 



(14) 



The definition of the running velocity defines the running coupling constant: 



ff*(k) 



47rufl(k) 



(15) 



The coupling constant appearing in the physical magnitudes is defined as a function of the Fermi velocity as g — t 
Two different renormalization prescriptions for the Fermi velocity measured at points k,A, ks are related by trie 
renormalization group equation: 



u(ki 



^lo K f> 



(16) 



This 



Hence, it is the RG what ensures that the exact theory is independent of the experimental point chosen in (|13[) . 
in turn guarantees the consistency of the renormalization procedure. 

The calculation in a dimensional regularization scheme follows exactly the same steps with an equivalent finite 
arbitrary part in the counterterm instead of F\ which is eliminated by the renormalization condition. That the 
procedure of renormalization does not depend on the cutoff is obvious considering that there exists a renormalization 
procedure (BPHZ scheme) that does not require the use any cutoff p8j . 

Notice that although being a tree level interaction, the counterterm is of order g. This diagram has to be added 
in the construction of the vertex functions at each given order in perturbation theory. In particular it will affect the 
second order diagrams computed in the next section. 

The one loop correction to the photon self-energy of Fig. [2](b), the polarization function, is finite and given by 



n(k,«) 



. e 



(17) 
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FIG. 4: Feynman diagrams contributing to the photon self-energy to second order in perturbation theory. 



This one loop result is independent on the nature of the interaction since only electron propagators appear in the 
calculation. Different interaction vertices describing Yukawa couplings, scalar potentials or disorder couplings may 
change the tensor structure but the diagram will remain finite. In QED(3+I) this diagram has a logarithmic singularity 
and the photon polarization has higher powers of logs at higher orders in perturbation theory what gives rise to the 
electric charge renormalization. 

Finally, the one loop vertex correction in Fig. [5] (c) is also finite and given by 

r° 1) (fc 0! k)~-| 7 e + finite. (18) 

The first term in this equation is the result of the explicit computation of the potentially divergent part of the diagram 
that turns out to be finite. The rest is encoded in the "finite" part that is a complicated expression of all the variables 
including Feynman parameters not important for the present discussion. 

This finishes the renormalization of the theory to first order in perturbation theory (one loop level). Any observable 
quantity at this level can be computed with the Feynman diagrams of Figs. Q] and [3] and will be finite and independent 
of the regularization procedure. 

Notice that at this order in perturbation theory there are no Coulomb interaction corrections to the optical con- 
ductivity. We will discuss the renormalization of this function to the next order in perturbation theory in the next 
section. 



C. Photon propagator to second order in perturbation theory 

In most of the standard perturbative renormalizable theory like QED(3+I), second order (two loops) corrections 
only change slightly the values of the observables computed at first order and do not change the physics set at 
first order. A different case occurs when a primitively divergent diagram appears at second order, the classical 
example being the two loops self-energy correction in A</> 4 in (3+1) dimensions that gives rise to the wave function 
renormalization absent at the one loop level [34| . In the graphene case the new physical feature that appears at 
second order is related to the electron wave function renormalization whose absence at the one loop level is due to 
the instantaneous nature of the photon propagator. Since the main issue of the present work is the determination of 
observables in the renormalized theory we will focus on the photon propagator that is the most important object in 
the transport properties. 

The full propagator at order g 2 will be given by the sum of the diagrams shown in Fig. |4j Diagrams Fig. |4]Ja) 
and Fig. [He) are finite in this model. The diagram in Fig. HJb) contains a sub divergent diagram which is identified 
as the first order correction to the electron self-energy diagram discussed in the previous section. This divergence is 
canceled by the diagram in Fig. Hfd) leaving a finite, unambiguous, result. 

The sum of diagrams 0|b) and0|d) reads: 



IT + rL 



d 3 k 

7>)3 



7 °G(fc) 



.(i) 



d 3 p 
J2nj~ 3 



7 G(fc+p) 7 ° 1 



G(fc) 7 °G(fc + (Z ), 



(19) 



where T^ t ' is given by eq. ©• We can identify the piece in the brackets as the renormalized electron self-energy 
with both contributions from the diagram in Figj2ja) and its counterterm shown in Fig|3] Hence the inner loop has 
already been computed, and the potentially divergent part that remains is: 



n 6 +n d 



32irv F J (2tt) 



d 3 k 7° (7% + v Fl • k) 7 • k (7°/c + vfI ■ k) 7° (7°(fco + q ) + «f7 ■ (k + q)) log £r 
' rTr , ,„ „ ^ ^ ^ (20) 



{-kl + V F k 2 n-(k + q,Y + V F (k + q)2) 
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which an explicit calculation shows to be finite. Thus, the renormalization procedure for the photon propagator is 
complete to second order in perturbation theory by means of the renormalization condition that made the electron 
propagator finite to first order in perturbation theory. 

The photon propagator is directly related to the optical properties of the system since it is included inside the 
definition of the dielectric function e(u, q) as: 

-^- = l + F (q)n(q, W ), (21) 

which in the limit u — ¥ gives the static screening properties of the system and is purely real. In undoped graphenc 
no finite screening length is generated and the effective Coulomb potential remains long ranged. Both the real and 
imaginary parts of the photon self-energy are related to physical observables. The imaginary part of the photon 
self-energy is related to the frequency-dependent conductivity through: 

Recr(w) = lim -^Imn(w,q). (22) 

q->0 q^ 

From which the optical response of the system can also be obtained. The non-interacting optical conductivity can be 
computed from the one-loop diagram in fig EJa) . This result is well known: 

2 

*o(w) = ~- (23) 

Therefore the optical conductivity at first order is a constant independent of the interaction. The Coulomb interaction 
corrections to the optical conductivity are obtained from eq. (|22[) plugging in the photon propagator at the two loops 
level. Since we have not computed the finite parts of the diagrams involved we can not give a precise number but from 
the discussion of the previous section it is clear that the result does not depend on the regularization prescription and 
it is determined by the chosen experimental data fixing the Fermi velocity. 



IV. CONCLUSIONS AND DISCUSSION 



In this work we have addressed the determination of physical observables in graphene in the presence of Coulomb 
interactions of density-density type modeled with a static Coulomb potential within a QFT perturbative renormaliza- 
tion scheme. We have shown that all the divergences up to two loops of the physical quantities are associated to the 
one loop electron self-energy and can be regulated by a proper renormalization of the Fermi velocity. The consistency 
of the scheme presented was exemplified by the renormalization of the photon polarization at the two loops level. We 
have shown that the Coulomb interaction corrections to the optical conductivity can be fixed unambiguously with a 
proper renormalization condition (experimental input) and do not depend on the regularization scheme. 

The QFT renormalization procedure outlined in this work shows that it can also be applied to condensed matter 
theories where there is an an ultraviolet physical cutoff (the inverse of the lattice spacing A ~ 1 /a) that prevents the 
appearance of infinities. Although initially devised to get rid of ultraviolet divergences, the QFT renormalization in 
more modern approaches is a way to define the physical parameters of effective theories in a given range of energies 
that do not depend on the details or parameters at higher energies. 

As we have seen the most important parameter in the computation of the observable quantities in graphene is the 
renormalized Fermi velocity that defines the Coulomb coupling constant. In the model discussed in this work the 
Fermi velocity grows without bound in the infrared (l2j ] . This sets a lower bound on the validity of the model that 
breaks down when the Fermi velocity approaches the speed of light c. The scale defines an infrared cutoff S for the 
theory, and can be computed as: 



fcfjexp 



167r(c - v(k R )) 



(24) 



This estimate depends on the renormalization point kn and on v(k R ). Plugging in the measured values we realize that 
the huge exponential suppression gives an infrared cutoff below any experimental resolution. This does not impose 
any real bound on the validity of the model from a experimental point of view. The limits of validity of the static 
model are set by internal consistency of the theory: the choice of a charge-charge interaction made to model Coulomb 
interactions is consistent with the static approximation since both are related to the ratio of Fermi velocity over 
the speed of light. When the Fermi velocity increases it would be more consistent to consider a retarded Coulomb 
interaction and a full interaction vertex. The retarded model was analyzed up to one loop in [l2j and has been revised 
recently in (35j | . 
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For the renormalization procedure to be complete, the fine structure constant of graphene ag = 47r e £t)F must be 
fixed, being and important quantity that appears in all the electronic properties of the system. A precise experimental 
determination is a necessity as the experiments in graphene are reaching a high degree of accuracy. The program to fix 
this constant can be similar to the one followed in the case of the electromagnetic fine structure constant [29| . In the 
case of graphene the theoretical determination is much simpler since there are no mass parameters relations involved 
in the calculations. Moreover it is very interesting the additional fact that the quantity determining uq (Fermi 
velocity) is by itself an observable related to the one particle properties of the system. We expect that a proper 
combination of photoemission .5, 36* optical [8j and transport [7|, |37J measures with the corresponding calculations in 
the renormalization scheme described here should be enough to determine olq as precisely as needed both theoretically 
and experimentally as it happens in QED(3+1). 

The renormalization program described in this work can be carried out to all orders in perturbation theory. As 
described in [12, [H| j the electron self-energy at the two loops level has a logarithmic singularity that induces a wave 
fuction renormalization. A new counterterm of order g 2 has to be added that affects the photon polarization function 
at the three loops level. This counterterm can be fixed by the same renormalization condition used to fix the Fermi 
velocity: requiring that the two point function at the two loops order has the same form as the free one, with a 
finite residue at kn given by a second experimental input. In the present model all physical quantities are fixed by 
renormalizing the electron propagator only. It is worth noticing that with the instantaneous Coulomb interaction 
discussed in this work, the wave function renormalization does not give rise to an anomalous exponent and the residue 
of the quasiparticle in the random phase approximation done in |14J is finite at the Fermi surface implying that, 
despite the anomalous lifetime [13j the system is closer to a Fermi liquid than to a marginal Fermi liquid [38] . 
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